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Laboratory studies reveal a deformational instability that leads to a drained region (dry spot) in an initially flat liquid layer 
(with a free upper surface) heated uniformly from below. This long-wavelength instability supplants hexagonal convection cells 
as the primary instability in viscous liquid layers that are sufficiently thin or are in microgravity. The instability occurs at 
a temperature gradient 34% smaller than predicted by linear stability theory. Numerical simulations show a drained region 
qualitatively similar to that seen in the experiment. 

PACS numbers: 47.20.Dr, 47.20.Ky, 47.54.-fr, 68.15.-be 


Benard’s observation in 1900 Q of hexagonal convec¬ 
tion patterns launched the modern study of convection, 
pattern formation, and instabilities; yet understanding of 
the surface-tension-driven regime in which Benard per¬ 
formed his experiments is still far from complete. Block 
and Pearson first showed how temperature-induced 
surface tension gradients (thermocapillarity) caused the 
instability observed in Benard’s experiments. Pearson’s 
linear stability analysis with a nondeformable liquid-gas 
interface yielded an instability at wavenumber q = 1.99 
(scaled by the mean liquid depth d) and Marangoni num¬ 
ber Me = 80, where M = ap ^ Td/pun (see Fig. ||) 
expresses the competition between the destabilization by 
thermocapillarity and the stabilization by diffusion. 

A deformable free surface allows a second type of 
primary instability in which a perturbation creates a 
nonuniform liquid depth and temperature profile [^,^. 
Thermocapillarity causes cool, elevated regions to pull 
liquid from warm, depressed regions (see Fig. |^). The 
instability appears with a long wavelength since sur¬ 
face tension stabilizes short wavelengths. Linear stability 
analyses that allow for deformation reveal this instability 
(see Fig. ||) at zero wavenumber {q — 0) and Me = 
[^,||, where the Galileo number, G = gd^/vK {g is the 
acceleration of gravity), gives the relative strengths of 
the stabilizing mechanisms of diffusion and gravity and 
thus determines which instability will form. For suffi¬ 
ciently thin, viscous liquid layers or small g (e.g., micro¬ 
gravity), |G < 80, so the long-wavelength mode should 
become the primary instability (see Fig. ^). This long- 
wavelength instability has not been experimentally inves¬ 
tigated. 

In this Letter we describe experimental observations 
of the onset of the long-wavelength instability and com¬ 
pare these observations to linear stability theory. The 
instability leads to a large-scale drained region with di¬ 
ameter ~ lOOd. The qualitative features of the instabil¬ 
ity are compared to nonlinear theory based on a long- 


wavelength evolution equation. We also explore the com¬ 
petition between the long-wavelength and hexagonal in¬ 
stabilities and study the physical mechanism that selects 
which pattern will appear. For a range of liquid depths, 
both patterns coexist. 
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FIG. 1. A surface-tension-driven Benard (Marangoni) con¬ 
vection cell (not to scale) contains both a liquid (silicone oil) 
layer of mean depth d and local depth h{x, y, t) and a gas 
(air) layer of mean thickness dg. The mean temperature 
drop across the liquid layer is AT. The liquid has density 
p (0.94 g/cm®), kinematic viscosity v, thermal dlffusivity k 
(0.0010 cm^/s), and temperature coefficient of surface ten¬ 
sion ctt = I da/dT \ (0.069 dynes/cm°C). 



FIG. 2. Marginal stability curves for layers of 10 cS silicone 

oil at two thicknesses: d = 0.1 cm (-), hexagonal convection 

cells form at Me = 80 with q = 1.99 (G = 10“*); d = 0.01 

cm (-), the long-wavelength {q = 0) instability forms at 

Me = 6.7 (G=10). 

We study a thin layer of silicone oil that lies on a 
heated, gold-plated aluminum mirror and is bounded 
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above by an air layer (see Fig. |^. A single-crystal sap¬ 
phire window (0.3-cm-thick) above the air is cooled by 
a temperature-controlled chloroform bath. The temper¬ 
ature drop across the liquid layer is calculated assuming 
conductive heat transport Q and is typically 0.5-5 °C. 
We use a polydimethylsiloxane silicone oil with a vis¬ 
cosity of 10.2 cS at 50 °C. The circular cell (3.81 cm inner 
diameter) has aluminum sidewalls whose upper surface is 
made non-wetting with a coating of Scotchgard. The ex¬ 
periments are performed with 0.005 cm < d < 0.025 cm 
and 0.023 < dg < 0.080 cm (typically dg = 0.035 cm); the 
corresponding fundamental wavevectors are in the range 
0.008 < q = 2Trd/L < 0.040 (^ 1). The liquid layers are 
sufficiently thin that buoyancy is negligible ^ . The gap 
between the sapphire window and the mirror is uniform 
to 1%, as verified interferometrically. The liquid surface 
is initially flat and parallel to the mirror to 1% in the 
central 90% of the cell, with a boundary region near the 
sidewalls due to contact line pinning. This initial depth 
variation is accentuated by thermocapillarity as AT is 
increased; measurements of depth variation show a 10% 
surface deformation (in the central 90% of the cell) at 
3% below onset. For visualization, we use interferome¬ 
try, shadowgraph and infrared imaging (256 x 256 InSb 
staring array, sensitive in the range 3-5 /im). 



FIG. 3. Dry spot in both experiment and numerical simu¬ 
lation. The experimental picture shows the measured bright¬ 
ness temperature as a function of position along the interface 
just above (~ 10%) onset; d = 0.011 cm. Numerical simula¬ 
tion of a long-wavelength evolution equation shows the depth 
as a function of position for 1% above onset of linear insta¬ 
bility. The depth of the liquid goes to zero in the drained 
region. 

Above a critical AT, the liquid layer becomes unsta¬ 
ble to a long-wavelength draining mode that eventually 
forms a dry spot [see Fig. ||(a)]. The drained region 
takes several hours (of order a horizontal diffusion time, 

Ik) to form. The dry spot is not completely dry since 
an adsorbed layer ~ 1 ^m remains fl^ . The size of the 
drained region is typically one quarter to one third the 
diameter of the entire cell. 

Our measurements for the onset of instability are com¬ 
pared to the prediction of linear stability theory §00 
in Fig. The results are given in terms of the dynamic 


Bond number, B = G/M = pgd^/axAT, which is the 
relevant control parameter for the long-wavelength insta¬ 
bility; T is a measure of the balance between gravity’s 
stabilizing influence and thermocapillarity’s destabilizing 
effect. Figure ^ shows that the measured critical values 
of B~^ are independent of liquid depth and viscosity, as 
predicted by theory. However, B~^ in the experiment 
is 34% smaller than predicted. We do not believe this 
discrepancy is due to systematic errors in the character¬ 
ization of our experiment {e.g., geometry, fluid proper¬ 
ties) since experiments in the same convection cell using 
thicker liquid layers find onset of hexagons in agreement 
with another experiment 111 and linear stability theory 

i- 



FIG. 4. Measurements of onset compared to the theoretical 

prediction, B~^ = | (-). (Solid circles, 10.2 cS silicone oil; 

solid triangles, 30 cS silicone oil.) The weighted mean of the 
data is — 0.44 ± 0.03 (-). There is an additional un¬ 

certainty of ±0.06 in Bc^ from the uncertainty in the thermal 
properties of the silicone oil. 

The long-wavelength and the hexagonal modes become 
simultaneously unstable at a critical liquid depth dc- 
Near this critical depth, both modes of instability com¬ 
pete and influence the formation of the pattern. The 
two modes are not mutually exclusive, but there is a fun¬ 
damental imbalance in their relationship: the presence 
of the hexagons suppresses the long-wavelength mode, 
while the presence of long-wavelength deformation may 
induce the formation of hexagons. Linear theory predicts 
dc = (120jzK/g)^/^ = 0.023 ± 0.001 cm; we observe the 
exchange of primary instabilities at dc = 0.025 ± 0.001 
cm. For d > dc this depth, hexagons are the primary in¬ 
stability [see Fig. |^(d)[; the hexagons smooth the large- 
scale temperature variations that would allow formation 
of the long-wavelength mode as a secondary instability 
when B~^ exceeds B~^. For d < dc, the long-wavelength 
mode is the primary instability [see Fig. 11(b)]. The fluid 
expelled from the forming dry spot increases the local 
height h{x,y,t) and thus the local Marangoni number 
[oc h(x, y, t)] in the newly formed elevated region; for 
0.017 cm < d < dc, hexagons form in the elevated re¬ 
gion since the local M in this region exceeds Me (80) for 
the onset of hexagons [see Fig. 11(c)] . For d < 0.017 cm. 
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hexagons do not form at the onset of the long-wavelength 
mode, but can form in the elevated region for AT suffi¬ 
ciently above onset. Similar mode competition phenom¬ 
ena have been studied theoretically for solutocapillary 
convection |^. When d dc, increasing AT above 
ATc increases the area of the dry spot; no qualitatively 
new structures are observed. At fixed AT, the dry spot 
is stable. 
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FIG. 5. Infrared images from experiments with increasing 
fluid depth (increasing G). (a) 10% below onset of instability 
(d — 0.033 cm), (b) For d — 0.011 cm, the long-wavelength 
mode is the primary instability. The dark region is the dry 
spot, (c) At d = 0.022 cm, the long-wavelength and hexago¬ 
nal modes coexist. A droplet (white circle) is trapped within 
the dry spot (dark oval); the liquid layer is strongly deformed 
(white annulus) between the dry spot and the hexagonal pat¬ 
tern. (d) For d > 0.025 cm, hexagons are the primary insta¬ 
bility (here, d = 0.033 cm). 

To understand the long-wavelength instability more 
fully, we study a long-wavelength evolution equation de¬ 
rived by Davis for an insulating upper boundary 

and later considered by other authors 00 for more 
general upper boundary conditions. The evolution equa¬ 
tion for the local liquid depth h{x, y, t) (scaled by d) with 
an insulating upper boundary is 

+V v^l =0 

where the domain of both x and y is [0,27r), time is scaled 
by d? /B is the dynamic Bond number as before, and 
the static Bond number is Bo = pgV^/a. The first term 
in curly brackets describes the effect of thermocapillar¬ 
ity; the second, gravity; and the third, surface tension. 
A linear stability analysis of this equation agrees with 


the linear stability analysis of the full fluid equations ref¬ 
erenced earlier. Numerically integrating the evolution 
equation in time using a pseudo-spectral code with pe¬ 
riodic boundary conditions (which automatically satisfy 
the incompressibility condition), we find that when the 
system is unstable, an initially infinitesimal perturbation 
grows until it forms a dry spot, as shown in Fig. |(b) {Bo 
= 1200 in the simulation, 700 in the experiment). Dry 
spots in the simulation do not saturate (unlike the ob¬ 
served spots), since the assumptions leading to the evolu¬ 
tion equation are violated once the dry spot has formed; 
in this regime, the simulation quickly loses spectral con¬ 
vergence and breaks down. The shape of the dry spot 
near onset is independent of parameter values except for 
To, which governs the size of the drained region. For 
large To, sharp structures can form; for small To, round¬ 
ing by surface tension prevents the formation of sharp 
structures. 

We find no stable, intermediate solution of the evo¬ 
lution equation between no deformation and a drained 
region. Our weakly nonlinear analysis of the evolution 
equation predicts a subcritical instability for all parame¬ 
ter values. We have performed one- and two-dimensional 
simulations as well as an analysis of the evolution equa¬ 
tion to look for any turnover in the backwards unstable 
(subcritical) branch of the bifurcation. For an insulating 
upper boundary, we have shown analytically that no sta¬ 
ble deformed solutions exist, in agreement with previous 
authors Q. The unstable branch continues backwards 
until T“~= i, at which point it disappears since the 
unstable solution past that point is unphysical (h < 0); 
nowhere does the backwards unstable curve turn over in 
a saddle-node bifurcation. Simulations using more gen¬ 
eral upper thermal boundary conditions pT| ] also find no 
stable deformed solutions. 

In conclusion, we have observed the formation of a 
large-scale dry spot due to a long-wavelength deforma- 
tional instability. The formation of such dry spots could 
be a serious problem in the planned use of liquid layers 
(even as thick as one centimeter) in microgravity environ¬ 
ments, where fluid motion is driven primarily by surface 
tension gradients (buoyancy effects are negligible). 

The structure of the dry spot agrees qualitatively with 
numerical simulations. Linear stability theory correctly 
predicts the functional dependence of the onset on ex¬ 
perimental parameters, but the onset in the experiment 
occurs at B~^ (or ATc) 34% smaller than predicted by 
linear theory. This discrepancy in the onset may be due 
to the difference in lateral boundaries between theory 
and experiment. For example, all theory to date has as¬ 
sumed periodic boundary conditions, which become sus¬ 
pect when the wavelength of the mode under consider¬ 
ation is of order the system size. The long-wavelength 
evolution equation is not valid near a boundary and it is 
not clear what boundary conditions other than periodic 
might be used for the evolution equation. Future nu- 
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merical simulations employing realistic, finite horizontal 
boundaries may answer this question. 

We thank S. H. Davis, R. E. Kelly, and E. L. 
Koschmieder for useful discussions. This research is sup¬ 
ported by the NASA Microgravity Science and Appli¬ 
cations Division (Grant No. NAG3-1382). S.J.V.H. is 
supported by the NASA Graduate Student Researchers 
Program. 


Swift and H. L. Swinney, submitted to Phys. Rev. Lett. 

[14] A. A. Golovin, A. A. Nepomnyashchy and L. M. Pismen, 
Phys. Fluids 6, 34 (1994). 

[15] S. H. Davis, in Waves on Fluid Interfaces, edited by R. 
E. Meyer (Academic Press, New York, NY, 1983), p. 291. 

[16] S. H. Davis, Annu. Rev. Fluid Mech. 19, 403 (1987). 

[17] B. K. Kopbosynov and V. V. Pukhnachev, Fluid Mech. 
Sov. Res. 15, 95 (1986). 

[18] A. Oron and P. Rosenau, J. Phys. II (France) 2, 131 
(1992). 


* Electronic mail: svanhook@chaos.ph.utexas.edu 

^ Electronic mail: schatz@chaos.ph.utexas.edu 

* Electronic mail: swinney@chaos.ph.utexas.edu 

[1] H. Benard, Rev. Gen. Sci. Pure Appl. 11, 1261 (1900). 

[2] M. J. Block, Nature 178, 650 (1956). 

[3] J. R. A. Pearson, J. Eluid Mech. 4, 489 (1958). 

[4] L. E. Scriven and G. V. Sternling, J. Eluid Mech. 19, 321 
(1964). 

[5] K. A. Smith, J. Eluid Mech. 24, 401 (1966). 

[6] M. Takashima, J. Phys. Soc. Japan 50, 2745 (1981). D. 
A. Goussis and R. E. Kelly, Int. J. Heat Mass Trans. 
33, 2237 (1990). G. Perez-Garcia and G. Carneiro, Phys. 
Fluids A 3, 292 (1991). 

[7] E. L. Koschmieder and M. I. Biggerstaff, J. Eluid Mech. 
167, 49 (1986). 

[8] We use a silicone oil that has been distilled once to re¬ 
move low vapor pressure components, which can con¬ 
dense on the cool, upper plate [see M. E. Schatz and K. 
Howden, Exp. Fluids (to appear)]. Infrared images are 
made using a mixture of 90% silicone oil and 10% 20 cS 
(at 50 °G) infrared-absorbing (4.61 /rm) polymethylhy- 
drosiloxane silicone oil. 

[9] Surface tension is the dominant driving force when 

M/R > 1, where the Rayleigh number is i? = 

agATd^/i/K and a is the volume coefficient of expansion. 
[See S. H. Davis and G. M. Homsy, J. Fluid Mech. 98, 
527 (1980).] Buoyancy is negligible in our experiments 
since M/R > 110. 

[10] J. P. Burelbach, S. G. Bankoff and S. H. Davis, Phys. 
Fluids A 2, 322 (1990). 

[11] This assumes an insulating upper boundary. The more 
general case is = |(1 + B), where the Biot number 
H characterizes the thermal properties of the liquid-gas 
interface. In the long-wavelength limit, H = kgd/kjdg, 
where k and kg are the thermal conductivities of the liq¬ 
uid and gas, respectively [see L. Hadji, J. Safar and M. 
Schell, J. Non-Equilib. Thermodyn. 16, 343 (1991)[. H 
ranges from 0.01 to 0.08 in our experiments. Values in 
Fig. ^ include the (1 -I- H) correction factor [that is, they 
are /{I + H)]. The pattern and the topology of the 
bifurcation diagram are independent of H. 

[12] These results are for an infinite system. A finite system 
has a lowest order mode of q — 2Tvd/L, which becomes 
unstable at a slightly higher (6%) value of the control 
parameter. We do not include this correction in Fig. ^ 

[13] M. F. Schatz, S. J. VanHook, W. D. McCormick, J. B. 


4 



